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Light-induced atomic desorption and diffusion of Rb from porous alumina 
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We present the first study of light induced atom desorption (LIAD) of an alkali atom (Rb) in 
porous alumina. We observe the variation due to LIAD of the rubidium density in a vapor cell 
as a function of illumination time, intensity and wavelength. The simple and regular structure of 
the alumina pores allows a description of the atomic diffusion in the porous medium in which the 
diffusion constant only depends on the known pore geometry and the atomic sticking time to the 
pore wall. A simple one-dimensional theoretical model is presented which reproduces the essential 
features of the observed signals. Fitting of the model to the experimental data gives access to 
the diffusion constant and consequently the atom-wall sticking time and its dependence on light 
intensity and wavelength. The non-monotonic dependence of the LIAD yield on the illumination 
light frequency is indicative of the existence of Rb clusters in the porous medium. 

PACS numbers: 68.43.Tj 66.30.Pa 78.67.Rb 47.61.-k 78.67.Bf 



I. INTRODUCTION 



In a spectroscopic cell containing an alkali atom 
vapor, a substantial fraction of the atoms are adsorbed 
on the cell walls. At steady state, the gas density is in 
equilibrium with the adsorbed atomic fraction. In some 
cells, depending on the cell material or coating, when 
the cell is illuminated with moderate intensity (1 - 1000 
mW/cm 2 ) nonresonant light, a significant increase in 
the atomic vapor density is produced as a consequence 
of the release of atoms from the cell surface into the gas 
phase. Such effect has been named light induced atomic 
desorption (LIAD) [l|. 

LIAD has received considerable attention in recent 
years due to its application as a light-controlled atom 
dispenser under high vacuum conditions. Such dispenser 
has been successfully used to load magneto-optical atom 
traps @, i, H| and hollow optical fibers @, S 0- Its 
use has also been considered for atomic magnetometers, 
gyroscopes and clocks @, H|. In addition, LIAD has 
attracted the attention of astrophysicists since it has 
been related to the observed abundance of alkaline 
elements in nonpermanent extraterrestrial atmospheres 

LIAD is understood as a non-thermal effect as opposed 
to light desorption produced with high power sources in 
which a significant heating results from light absorption 
by the substrate. In poly-dimetilsiloxane (PDMS) [ill ], 
paraffin 12 1 and sapphire [l3| a frequency threshold in 
the infrared, similar to that of the photoelectric effect on 
metals, has been observed. Also, an increasing efficiency 
of LIAD with light frequency has been reported in 
several samples [H, EH Ell ■ 

LIAD was first observed in sodium vapor glass cells in 
which the inner cell walls were coated with a thin layer 
of PDMS. The effect was also observed with K [lj|, 
Rb and Cs atoms [12j (sometimes in the presence of a 
buffer gas). Initially, it was considered that LIAD was 
specific to PDMS coatings [ll|. However, LIAD was 



later-on reported in cells coated with different polymers 
such as octadimethyl-cyclotetrasiloxane (OCT) [17| and 
paraffin [T2|. LIAD has also been observed on several 
uncoated surfaces such as glass [3, 0] , stainless steel Q 
and sapphire [l8| . More recently, LIAD has been studied 
in porous amorphous materials such as porous silica [lj| . 
All observations of LIAD in porous or coated surfaces 
present some common features such as the characteristic 
time scale of the atomic desorption (several seconds). 
However, other aspects may vary significantly between 
different atomic species and coatings depending also 
in the cell preparation procedure. In particular, large 
variations are observed in the desorption yield. In cells 
coated with PDMS, LIAD may result in an increase 
of the atomic gas phase density of several orders of 
magnitude [Tl| while only density-increase-factors of a 
few units were reported for paraffin 12j . Smaller factors 
were observed on uncoated surfaces as in the present 
study. The question on whether there is a common 
mechanism underlying all LIAD observations is still 
open [20| . 

The first tentative explanation of LIAD at the mi- 
croscopic level was suggested by Xu et al ll]. The 
mechanism involves the modification by light of the 
weakly bonded chemical complex formed between a 
PDMS molecule and the Na atom or Na2 molecule. 
More recently, this mechanism was further investigated 
through the measurement of the thermal distribution 
of desorbed atom velocities [2l|. This interpretation of 
LIAD is consistent with the observation of a threshold 
light frequency for LIAD in PDMS but leaves unex- 
plained several aspects of its dynamics. As discussed by 
Atutov et al [ljj], in addition to the atomic desorption 
from the surface, the diffusion of the atoms within the 
surface coating plays an essential role in the temporal 
evolution of LIAD. To a large extent, LIAD in coated 
surfaces is a consequence of light-induced modification 
of the atomic mobility and diffusion within the coating 
polymer. Atutov et al [l4[ have modelled such process 
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assuming a phenomenological dependence of the diffu- 
sion coefficient on light. Alexandrov et al have described 
the LIAD dynamics with the help of rate equations with 
a light dependent term representing the flux of atoms 
from the coating into the gas phase [T^. Recently, 
the model of LIAD in coated surfaces suggested by 
Atutov has been revised and improved by Rebilas and 
co-workers |22|, |23[ . 

LIAD from uncoated dielectric surfaces, such as sapphire 
or glass, deserves special consideration. In such systems, 
alkaline atoms can be individually adsorbed on the 
surface or agglomerated into clusters. The presence of 
clusters may result in a visible change of the sample 
transparency or even in coloration [24(. Blue-green 
coloration by Rb of otherwise transparent (or white) 
samples has been observed in several experiments 
including the ones described here. The role of the light 
in these samples is double since it can produce the direct 
desorption of the atoms from the dielectric surface and 
the evaporation of the atomic clusters [25|, [26|. Also, 
under suitable conditions, the light may also control 
the growth of the clusters from atoms in the vapor 
phase. Such conditions are favored in porous media 
where the desorbed atoms remain confined and available 
to participate in the cluster regrowth. A characteristic 
feature of the LIAD involving cluster evaporation, is the 
non-monotonic dependence of the desorption yield on 
the light frequency. Such behavior is interpreted as the 
consequence of resonant surface plasmon excitation in 
the cluster [25] . A second characteristic of these systems 
is its "memory" . The response strongly depends on the 
illumination history including the timing of the bright 
and dark periods and the corresponding color sequence 

A common feature governing the dynamics of both, 
the LIAD in polymer coated surfaces and in porous 
dielectrics, is the successive occurrence of two distinctive 
processes: i) atomic desorption from the surface (or 
cluster) ii) diffusion in the intermediate medium (either 
the polymer or the porous medium) prior to the atom 
release in the vapor phase. A precise modelling of LIAD 
should involve the simultaneous account of these two 
processes. The desorption mechanism is at present 
only qualitatively understood [ll], [2l|, [2f|. Also, little 
understanding is currently available on the mechanisms 
determining the variation with light of the atomic mo- 
bility in the polymer coating. Diffusion in porous silica 
is presumably simpler since the atomic motion inside 
the pores may be assumed to occur in a diluted vapor. 
However, the random nature of the pore geometry 
complicates the modelling of such process. 
In this paper we present the experimental study and 
theoretical modelling of LIAD with Rb atoms adsorbed 
in thin membranes of porous alumina. The porous 
medium is produced by anodization of aluminum and 
results in a very regular array of cylindrical pores with 
small size dispersion. The average diameter of the pore 
tube is 200 nm and its length 60 fim. In consequence, 




FIG. 1: (Color online) Experimental setup. Inset: Schematic 
cut of the alumina membrane for a plane parallel to the pores. 
BS: beam splitter, PD1, PD2: photodetectors. 



the pore geometry is well known and this allows a 
simple and accurate description of the atomic diffusive 
motion in the pore. The diffusion process is the result 
of a random sequence of atomic free flight of the atoms 
confined within the pore wall. After a collision with the 
wall, the atom sticks to the wall for some time after 
which it is desorbed again flying with random direction 
and velocity. Under such assumptions, the diffusive part 
(ii) of the LIAD dynamics, including the release of atoms 
into the gas phase, is determined by the pore geometry 
and the atom-surface sticking time. Such picture allows 
us to model the atomic diffusion in the pores as well 
as the atom exchange between the porous material and 
the outside gas volume. The resulting LIAD dynamics 
is mainly dependent on the atomic diffusion constant 
in the porous medium which can be obtained from the 
fitting of the observed variation of the atomic density 
in the gas phase. From the knowledge of the diffusion 
constant, the atom- wall sticking time can be determined. 
We have investigated the variation of the sticking time 
with light intensity and color. 



II. EXPERIMENT 

We have used porous alumina membranes manufac- 
tured by Whatman International Limited. The circular 
flat membranes have a diameter of one inch and a thick- 
ness of 60 /Ltm. The membrane is traversed by a regular 
array of cylindrical hollow tubes with 200 nm diameter. 
The tubes form a honey comb like array with a pore den- 
sity of 10 9 /cm 2 . The diameter of the pores are uniform 
over most of their length. On one side of the membrane, 
along I /im, the pores divide into several smaller branches 
with 20 nm typical diameter (see the inset in Fig. [TJ. 
Before contact with the Rb vapor the porous membranes 
are translucid and white. In order to fit into the vacuum 
glass cell, the membranes are divided in pieces of typi- 
cally 0.5 cm 2 . 

The experimental setup is shown in Fig. [T] We have 
used a vacuum glass cell with 2.5 x 2.5 x 4.5 cm dimen- 
sions. The cell is connected via a glass-to-metal tran- 
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sition fitting to an ion pump and a metallic Rb reser- 
voir. The Rb density in the glass cell is monitored by 
measuring the absorption of a laser beam issued from an 
extended cavity diode laser. Using a saturated absorp- 
tion setup, the laser frequency is stabilized to the 85 Rb 
F = 3 i-> F' = 4 transition in the D2 line (780 nm). In 
order to increase the absorption signal, the laser beam 
crosses the cell several times. We have used a balanced 
detection scheme to reduce sensitivity to laser intensity 
fluctuations. Half of the laser power is sampled before 
the cell and detected with a photodiode. A second pho- 
todiode monitors the intensity of the beam transmitted 
through the cell. The outputs of the two photodiodcs 
are subtracted. In order to eliminate noise from ambi- 
ent light, including the light used for the LI AD, the laser 
beam is modulated with a chopper and lock-in detected. 
The illumination of the porous alumina samples is made 
with high power LEDs (100 mW) in order to have a non- 
thermal source with a well defined spectrum. Three dif- 
ferent LEDs were used centered at 455, 504 and 617 nm 
(typical spectral width 10 nm). An optical arrangement 
(not shown in Fig. [1} allows a uniform illumination of 
the porous sample by the LED light. 



Prior to the introduction of the porous alumina 
membrane, the glass cell was evacuated (10 -6 torr) and 
baked for several hours at 300 C. Such precaution ap- 
peared to be essential since we have observed significant 
LIAD from the unbaked cell presumably due to some 
uncontrolled coating. After the baking procedure, the 
LIAD from the cell walls was negligible. Following the 
cell cleanup, several pieces of the alumina membrane 
were introduced and vacuum baked for several days 
at 150 Celsius. The pieces of alumina lied on the cell 
bottom. We had no control on the side of the membrane 
that faces the cell wall, so some of the pieces present 
the largest pore apertures toward the cell bulk volume 
while others present the narrow ramification ends. After 
the initial cleanup of the alumina, the cell was returned 
to room temperature and the valve separating the cell 
from the metallic Rb reservoir opened. Keeping the 
Rb reservoir and the vacuum connecting tubes slightly 
heated (~ 50 C), the Rb was allowed to diffuse into the 
cell and the porous alumina. After a few days, a visible 
blue coloration appeared in the alumina indicating the 
presence of Rb. After a sufficiently long period all the 
samples were dark blue. However two different blue 
tones were observed among the samples. We interpret 
such difference as a consequence of the two possible 
orientations of the membrane pieces with respect to the 
cell wall. The blue coloration is an indication of the 
formation of Rb clusters [24J. We have checked that 
the cluster formation is entirely reversible. The original 
white coloration of the alumina could be recovered after 
pumping the cell during a few hours while illuminating 
with an incandescent lamp. 
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FIG. 2: Typical observed variations of the relative atomic 
vapor density for two different illumination-time intervals. 



III. EXPERIMENTAL RESULTS 

We have observed the LIAD of Rb from the porous alu- 
mina by monitoring the laser absorption in the cell bulk 
while turning on and off the illumination by a LED. We 
have recorded the relative variation of the vapor density 
S(t) = (p(t) — Po)/po as a function of time where p(t) 
is the density of Rb in the cell and po the equilibrium 
density in the dark. Fig. [5] shows two typical records ob- 
tained with the same illumination for two different light- 
on intervals (500 and 600 s). In general, the relative den- 
sity reaches a maximum S max after a few tens of seconds 
depending on light intensity. After that the Rb density 
slowly decreases towards a new steady state in the pres- 
ence of light. When the light is turned off, 5 decreases 
on a time scale comparable to the rise time. Two differ- 
ent behaviors have been observed for long times after the 
light switching off. Either the density returns monotoni- 
cally to the initial equilibrium density po or drops below 
Po by an amount e (see Fig. [5]) after what it slowly grows 
towards po . The later behavior is observed if the light in- 
tensity and the illumination interval are sufficiently large. 

Fig. [3] shows the evolution of the Rb density in the 
cell for two different illumination intensities. Notice 
the variation in the shape of the trace. Similar shape 
variations were also observed in experiments with porous 
silica 19], although not reproduced by the proposed 
theoretical models. 

We have observed that the efficiency of the LIAD pro- 
cess depends on the porous alumina history, as was also 
noticed in other systems [l2|, [H| . A monotonic reduction 
in the maximum relative Rb density variation S max is ob- 
served for several successive illumination cycles keeping 
constant the light intensity. In addition, as the inten- 
sity is changed between successive illumination periods, 
the signal variation is different depending on whether the 
light intensity is increased or decreased (see Fig. Q|. For 
low enough light intensities the system is not appreciably 
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FIG. 3: (Color online) Two records of the relative atomic 
density variation illustrating the difference in shape for dif- 
ferent illumination intensities. The traces have been re-scaled 
to signal the difference in shape. 
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FIG. 4: Maximum relative atomic density S ma x as a function 
of the illumination intensity. The measurements were regis- 
tered with a sequence of illumination intervals of 140s followed 
by intervals of 600s without illumination. Solid (hollow) tri- 
angles: increasing(decreasing) illumination intensity. 



modified by the illumination and a linear dependence of 
the LIAD yield on light intensity is observed. The non- 
linear dependence, visible in Figf?] for large intensities, 
can be attributed to the depletion of the available Rb 
inside the nano-pores. 
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FIG. 5: (Color online) Maximum relative atomic density 5 ma x 
as a function of the illumination intensity for different illumi- 
nation wavelengths (455, 505 and 617 nm). 



wavelengths. In FigOthe non- linear variation of S rnax is 
only noticeable for the highest intensities. 

From the linear fit of the data in Fig. [5] one can evalu- 
ate, for each wavelength, the coefficient a\ = 8 ma x 
proportional to the LIAD desorption rate per photon 
flux. We observe that this coefficient for blue, green 
and red light varies in proportion to 1, 0.73 and 1.1 
respectively. Such result indicate a non-monotonic 
evolution of the LIAD yield with photon energy. 



IV. THEORETICAL MODEL 

We model the evolution of the atomic density inside the 
cylindrical pores as a one-dimensional diffusion process. 
The typical sticking time t s of alkali atoms on dielectric 
surfaces is of the order of tens to hundreds microseconds. 
After desorption, the atoms leave the internal surface of 
the pore with thermal velocity in a random direction with 
a Lambertian probability distribution [27| . The gas den- 
sity inside the pores is considered sufficiently low to ne- 
glect the collisions between flying atoms. At room tem- 
perature and for hundred nanometers tube diameters, af- 
ter a few nanoseconds flight, the atom is again adsorbed 
on the pore surface. Since the pore length is much larger 
than its diameter, we can consider that the atoms exe- 
cute a one dimensional random walk, along the pore axis, 
characterized by the diffusion constant (see Appendix) : 



I 2 d 2 
D = — = — 
2r 3r 



(1) 



Fig. [5] shows 5 max as a function of illumination 
intensity for three different wavelengths. The measure- 
ments were taken alternating the three available light 
colors successively for each intensity. The effect of the 
history on the LIAD efficiency is so reduced for the 
comparison among measurements taken with different 



Where I 2 is the mean square displacement per step in the 
random walk, r is the mean interval between steps which 
is essentially determined by the sticking time r ~ r s on 
the internal pore surface and d is the pore diameter. 

The atomic desorption is described by a reduction of 
t s induced by the light. We assume a simple linear de- 
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FIG. 6: a) Schematic one-dimensional model for the system, 
b) Illustration of the atomic release form the pore end into 
the gas cell (J + in Eq. [8$. 



pendence: 



T s n(l - kI) 



(2) 



where t s q is the sticking time in the dark, / is the light 
intensity and k a coefficient which is wavelength depen- 
dent. In consequence: 



D = 



Do 



(1 - kI) 



(3) 



with D Q = d 2 /(3T s0 ) being the atomic diffusion constant 
in the dark. 

Fig. [6] present a scheme of the one-dimensional model 
of the system. The cylindrical pore, considered closed 
on its left end has a total length L. To the right of the 
pore, the atomic vapor cell, associated to a length L , 
is connected to a reservoir accounting for the vacuum 
pumping and the external Rb reservoir. 

The (linear) density of atoms /x(y, t) inside the pore is 
described by the diffusion equation 



2t = D — 

dt dy 2 



(4) 



where y is the position coordinate inside the pore (see 
Fig®. 

The total number of atoms N in the cell is N — N g + 
N w where N g represents the atoms in the gas phase and 
N w the atoms adsorbed to the cell wall. The fraction of 
atoms in the gas phase relative to the total number of 
atoms is assumed to be a constant for given temperature 
and illumination conditions 281: 
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wavelength. The evolution of the atom number N in the 
cell is described by the equation 



dN 
~~dt 



dt 
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where J is the net atomic flux at the pore- vapor interface. 
The rate 7 describes the return to the equilibrium atom 
number N g o determined by the external pumping system 
and Rb reservoir. 

We separate the net flux J into two contributions J = 
J + + J~ describing the atoms leaving and entering the 
pore respectively. The flux of atoms entering the pores 
from the cell gas is given by: 



J- 



2L, 



(7) 



where v = (\v y \) is the mean magnitude of the atomic ve- 
locity in the direction of the pore. The simple geometry 
of our system allows the evaluation of J + without addi- 
tional assumptions by considering that the atoms within 
a mean step length I from the pore end have a probability 
1/2 for leaving the pore in the time interval r (see Fig. 
[6]b), then 
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The equations describing the evolution of the atomic den- 
sities fi and n = N g /L inside the pores and in the cell 
gas phase respectively are: 
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where we have introduced the parameters a = -, — ^7 — 

L c = (L + A) and 7 = ^^r- n is the equilibrium value 
of the atomic density in the vapor cell. 

The steady state densities in the dark inside the pores 
Ho and in the vapor cell no are linked through the condi- 
tion: 
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The boundary conditions at the pore ends are (see Eqs. 
Hand 



Here A represents an effective cell length corresponding 
to adsorbed atoms. Since the sticking time of the atoms 
to the cell walls can depend on light intensity, we consider 
that A depends on the illumination in the form: A = 
Arj(l — C-0 where £ is a coefficient that can depend on 
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From the above equations, it is possible to derive an 
approximate relation between the observed variation of 
the gas density in the cell and the corresponding change 
in the diffusion constant inside the pores. For this we no- 
tice that in our system, the return to equilibrium (under 
constant illumination) occurs on time scale which is long 
compared to the observation time. One can then con- 
sider that during the LIAD the total atomic population 
(inside the pores and in the cell) remains approximately 
constant: 

I v 

fiL + nL c ~ {——L + L c )n (12) 
ZL>o 

where we used Eq. 1101 

When the sample is illuminated, the LIAD effect re- 
distribute the atoms along the pore in a characteristic 
diffusion time L 2 /2D. If we assume that the gas phase 
density reaches its maximum n max in a time which is 
long compared to the diffusion time, one can consider 
that the corresponding atomic density inside the pores is 
approximately uniform [i ~ fj. m ax- Using Eqs. [5] (with 
7 = 0) and [T2l we have: 



an indication of cluster formation and cluster- light inter- 
action. The investigation of cluster formation is outside 
the scope of this work. 

We notice that the value of Lq in Table Q] is large com- 
pared to the length (< 10m) estimated from the actual 
glass cell volume. However, the total effective volume 
available to the atoms outside the porous alumina also 
depends on the vacuum system tubes and surfaces 0]. 
The parameter a reflects the dependence of the effective 
cell length on the illumination. The numerical fitting is 
quite insensitive to this parameter giving a large scat- 
tering of the results. The uncertainty in this parameter 
prevents the determination of a wavelength dependence. 
On the other hand, different values of the coefficient k 
are obtained depending on the wavelength of the illumi- 
nating light. 

From the parameters in Table [TJ we can check that the 
assumptions made for the derivation of Eq. fTS] are rea- 
sonable for our system. The estimate of the maximum 
relative vapor density variation obtained using Eq. [lj|] 
only differs in a few percents from the value resulting 
from the numerical integration of Eqs. [9] 
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Eq. [T2] can be used for a quick estimate of the relative 
variation of the diffusion constant from the observed 
change in the vapor density provided the second term 
inside the brackets in Eq. [13] is small. In the condi- 
tions of our experiment such term is of the order of unity. 

Some of the parameters appearing in the model can 
be directly determined for our system. From the porous 
alumina manufacturer we know that L = 60um, and 
d = 200nm. In consequence, (^ 2 ) = |gP m 1.6 x 10~ 7 m 2 . 
The mean atomic velocity at room temperature is 
v w 140m/s. The other parameters are determined 
through least square fitting of the numerical model to 
the experimental data. For this, we have numerically 
integrated the differential equations [5] with the boundary 
conditions given in Eqs. 1111 

Fig. [7] shows a typical experimental register together 
with the corresponding signal calculated from the model. 
The values of the parameters obtained from the fitting 
are presented in Table. [TJ The given uncertainties cor- 
respond to the scattering of the results of the fitting for 
different experimental runs. The value of Dq given in Ta- 
ble [IJ results from the average of the data obtained with 
all three excitation wavelength. Interestingly enough, the 
plot of the fitted values of Dq for different runs reveal a 
systematic grouping for each of the three colors used for 
LIAD (see Fig. [5}. In our model, Dq corresponding to 
the diffusion constant in the dark, is taken as constant 
and independent of the desorbing light color. However, 
the grouping observed in Fig. [5] may reveal a dependence 
of Dq on the illumination history. Such feature could be 
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FIG. 7: (Color online) Example of the fitting of the calculated 
signal (solid line) to the experimental data (circles). 

A comparison of the predictions of the theoretical 
model with the experimental observations is given in Fig. 
[9] for traces obtained with blue desorbing light. Except 
for the largest intensities, where the effects of saturation 
and illumination history are expected to be significant, 
the model correctly describes the growth of the LIAD 
signal with light intensity. Similar results arc obtained 
for the other colors used for illumination. 

Our theoretical model appears to correctly account for 
several features of the experimental signal. As shown in 
Fig. [10] the signal shape variation as a function of the il- 
lumination time interval is well described. In particular, 
the "undershoot" e of the vapor density below the initial 
density is well reproduced. Such "undershoot" is due to 
the small variation of the total number of atoms (due to 
the external pumping system) during illumination. As 
the illumination is turned off, the atoms are rapidly re- 
adsorbed by the porous alumina in a time shorter than 
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TABLE I: Fitted values of the parameters of the model. 
D [m 2 .s _1 ] L [m] li^ 1 ] crfmW^ 1 ] /c red [mW _:L ] «: 9ree „[mW _1 ] Kbiue [mW 1 ] 

2.8 ±0.5 x KT 11 106 ± 33 2.4±1.5xl0~ 4 2.2 ± 2.0 x 10~ 3 7.8 ± 0.7 x 10~ 3 4.6 ± 0.5 x 10~ 3 5.2±0.5xl0~ 3 
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FIG. 8: (Color online) Values of the diffusion constant Do 
obtained from the fitting of different experimental traces ob- 
tained with three illumination colors. 



the one required to equilibrate the cell with the exter- 
nal pump and Rb supply. The model also reproduces 
the difference in shape of the temporal evolution between 
low and large illumination intensities as shown in Fig II II 
Such shape variations were previously observed but not 
reproduced by existing models 

From the measured value of Dq using Eq. [1] one can 
determine the sticking time r s of the atoms to the pore 
walls. The obtained value r s o — 500/zs lies within the 
range of previous observations for alkali atoms on dielec- 
tric surfaces. A summary of the sticking times reported 
in the literature for several alkali atoms and surfaces is 
presented in Table [TTJ The value of t s is several orders 
of magnitude larger than the mean time-of-flight of the 
atoms between collisions with the pore walls To ~ Ins. 
At a given time, the fraction of atoms in the gas phase 
inside the pores relative to the total number of atoms 
participating in the diffusion is tq/t s ~ 10~ 5 . From the 
values in Table [J we estimate a relative variation of the 
atomic gas density within the pores of 60% for illumina- 
tion with 50 mW/cm 2 of red light. 



V. CONCLUSIONS 

We have studied LIAD of Rb atoms contained within 
alumina nanopores. We observed, as a function of time, 
the variations of the Rb density in the cell surrounding 
the porous alumina as illuminating light with different 
colors is turned on and off. We have shown that the 
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FIG. 9: (Color online) Comparison of the experimental data 
(circles and squares) with the calculated signal (solid lines) 
obtained with the parameters in Table U for A = 455nm 



observed signal evolution is determined by the diffusive 
motion of Rb atoms within the porous medium. Our ob- 
servations are consistent with the picture of atoms under- 
going a one dimensional random walk along the porous 
axis. Taking advantage of the well characterized geom- 
etry of the porous medium, a simple relation of the dif- 
fusion coefficient with the pore diameter and the atom- 
wall sticking time was established. Also, at the pores 
ends, the atom exchange between the gas cell and the 
porous medium is directly linked, without additional as- 
sumptions, to the parameters of the diffusive motion (Eq. 
U). 

The measurement of the diffusion constant gives di- 
rect access to the mean time between steps. This time 
is essentially a sticking time as the atoms remain most 
of the time absorbed to the pore wall. Our results indi- 
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FIG. 10: (Color online) Observed (circles and triangles) and 
calculated (solid line) signals for two illumination intervals. 



TABLE II: Reported alkali 
at room temperature 



dielectric surface sticking times 



Atom-Surface 


Sticking time 


Comments 


Cs-Pyrex 


1400^3 


1281 


Cs- sapphire 


< 160ms 


f28] 


Na-glass 


130ais 


[29] 


Rb- alumina 


500ais 


This work 



APPENDIX A: CALCULATION OF THE 
DIFFUSION COEFFICIENT IN A CYLINDRICAL 
PORE 

For a one dimension random walk in the direction y, 
assuming that the length and duration of the random 
step s are uncorrelated, the diffusion constant is given by 

m 




2t 2t 



(Al) 



where l y is the single step displacement in the direction 
y, v y is the y component of the particle velocity and t the 
time-of-flight of a given step, r is the mean time interval 
between steps. 




FIG. 11: (Color online) Calculated traces for different illumi- 
nation intensities. The plots have been re-scaled to signal the 
difference in shape. 



cate a linear decrease of the sticking time with the illu- 
minating light intensity for low light intensity. In addi- 
tion the sticking time modification appears to be depen- 
dent of the illuminating light frequency. The LIAD yield 
does not vary monotonically with light frequency for the 
three wavelength used. This suggests that the atom re- 
lease takes place, at least in part, from rubidium clusters 
where surface plasmon resonances contribute to the light 
absorption spectrum [24|, [25| . 



FIG. 12: (Color online) a) Cylinder and coordinate system 
considered in the calculation, b) Cross-section along the x, z 
plane. 

We consider particles free flying within the inner sur- 
face of cylinder with diameter d. A particle leaving the 
cylinder wall has a velocity given by: 



v z 

Vy 



v cos(0) 
vsm(6) sin(</>) 
v sin(0) cos(<j>) 



(A2) 



where v is the velocity modulus. See Fig. [T2] for angle 
definitions. 

The time-of-flight is given by: 
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V x ,z 



(A3) 



where l x g and v XiZ are the projections of the particle 
displacement and velocity over the x, z plane. We have: 
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lx,z — dcos(a) — d 



Using Eqs. [S2J [A3] and \A£ we obtain: 



dcos(9) 



v (cos 2 (6>) + sin 2 (6>) cos 2 ((/>)) 



and 



dsm(9) cos(9) sin(</>) 
(cos 2 (0) + sin 2 (0) cos 2 (</>)) 



(A4) 



(A5) 



(A6) 



The angular (Lambertian) distribution of the atoms 
leaving the surface is given by [13] : 



P(n)dfl = cos(d)d(l 



(A7) 



where is the solid angle. The thermal distribution for 
the magnitude of the atomic velocity is (27j : 



2 \k R T 



v exp 



(A8) 



with v r 



2k B T 



Using IA61 and [ATI after integration one gets: 



(K) = 3^ (A9) 
In a similar way, from IA5| IA7I and IA8I we obtain: 



T = (t)=d 



27TTO 



(A10) 



In our system, the time interval between flights is de- 
termined by the atom sticking time r s ( t ~ r s 3> tq). 
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